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Abstract 

Anisotropic flows (i?2 and V4) of light nuclear clusters are studied by Isospin-Dependent Quan- 
tum Molecular Dynamics model for the system of 86 Kr + 124 Sn at intermediate energy and large 
impact parameters. Number-of-nucleon scaling of the elliptic flow (1*2) are demonstrated for the 
light fragments up to A = 4, and the ratio of v±jv\ shows a constant value of 1/2. In addition, the 
momentum-space densities of different clusters are also surveyed as functions of transverse mo- 
mentum, in-plane transverse momentum and azimuth angle relative to the reaction plane. The 
results can be essentially described by momentum-space power law. All the above phenomena 
indicate that there exists a number-of-nucleon scaling for both anisotropic flow and momentum- 
space densities for light clusters, which can be understood by the coalescence mechanism in 
nucleonic degree of freedom for the cluster formation. 
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Anisotropic flow is of interesting subject 
in theoretical and experimental investi- 
gations on nuclear reaction dynamics in 
intermediate and high energy heavy-ion 
collisions [1,2,3,4,5,6,7,8,9,10]. Many stud- 
ies of the dependences of the 1-th and 2-nd 
anisotropic flows (the directed flow and el- 
liptic flow, respectively) on beam energies, 
mass number or quark number, isospin 
and impact parameter revealed much in- 
teresting physics about the properties and 
origin of the collective flow. In particular, 
recent ultra-relativistic Au + Au collision 
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experiments demonstrated the number of 
constituent-quark (NCQ) scaling from the 
transverse momentum dependent elliptic 
flow for the different mesons and baryons at 
the Relativistic Heavy Ion Collider (RHIC) 
in Brookhaven National Laboratory [11], 
it indicates that the partonic degree of the 
freedom plays a dominant role in forma- 
tion of the dense matter in the early stage 
of collisions. Several theoretical models 
have been successfully proposed to inter- 
pret the NCQ-scaling of hadrons at RHIC 
[12,13,14,15,16]. In these studies, a popular 
interpretation is assuming that the mesons 
and baryons are formed by the coales- 
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cence or recombination of the constituent 
quarks. Moving to the intermediate energy 
heavy ion collisions, the coalescence mech- 
anism has been also used to explain the 
formation of light particles and fragments 
[17,18,19,20,21]. In these studies, however, 
most observables which one focused on are 
the spectra of kinetic energy or momen- 
tum of light particles. While a few studies 
investigated the mass dependence of the 
directed flow [22,23]. Systematic theoret- 
ical studies on the flow and momentum 
space densities of different fragments in in- 
termediate energy domain in terms of the 
coalescence mechanism are still rare. In 
this context, we survey, for the first time 
to our knowledge, the nucleon number de- 
pendence of the anisotropic flows vi and 
t>4 in the intermediate energy heavy ion 
collisions with the coalescence scenario. 
Moreover, the power law behaviors of light 
fragment production in momentum space 
are also explored in our model simula- 
tion. Note that this power law behaviors 
have been experimentally demonstrated 
in violent collisions at the beam energies 
between 0.1 A and 15 A GeV [24]. In this 
Letter we use Isospin-Dependent Molec- 
ular Dynamics (IDQMD) model to simu- 
late 86 Kr + 124 Sn at 25 MeV/nucleon and 
larger impact parameters (b = 7 - 10 fm), 
and investigate the outcome of nucleonic 
coalescence mechanism on the anisotropic 
flows and momentum space densities of 
light particle production in intermediate 
energy heay ion collisions. 

Anisotropic flow is defined as the differ- 
ent ra-th harmonic coefficient v n of an az- 
imuthal Fourier expansion of the particle 
invariant distribution [2] 

dN 00 

— oc 1 + 2 V v n cos (n<f)), (1) 
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where is the azimuthal angle between the 
transverse momentum of the particle and 
the reaction plane. Note that in the coordi- 
nate system the z-axis along the beam axis, 
and the impact parameter axis is labelled 
as x-axis. The first harmonic coefficient v± 
represents the directed flow, v\ = (coscp) = 

(^), where p t = ^jp 2 x + p 2 is transverse 



momentum, v-i represents the elliptic flow 
which characterizes the eccentricity of the 
particle distribution in momentum space, 



V2 = (co S (20)) = (^-^), (2) 
Pt 

and t>4 represents the 4-th momentum 
anisotropy, 



pi 
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The model we are using in the present 
work is based on the Quantum Molecular 
Dynamics (QMD) approach which is an 
n-body theory to describe heavy ion re- 
actions from intermediate energy to 2 A 
GeV. It includes several important parts: 
the initialization of the target and the pro- 
jectile nucleons, the propagation of nucle- 
ons in the effective potential, the collisions 
between the nucleons, the Pauli blocking 
effect and the numerical tests. A general 
review about QMD model can be found in 
[25]. The IDQMD model is based on QMD 
model affiliating the isospin factors, which 
includes the mean field, two-body nucleon- 
nucleon (NN) collisions and Pauli blocking 
[26,27,28,29,30]. 

In the QMD model each nucleon is repre- 
sented by a Gaussian wave packet with a 
width \J~L (here L = 2.16 fm 2 ) centered 
around the mean position fi(t) and the 
mean momentum pl(t), 
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The nucleons interact via nuclear mean 
field and nucleon-nucleon collision. The 
nuclear mean field can be parameterized 
by 

U(p,r z )=a(^)+P(^y + l(l-T z )V c 
Po Po 2 
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with po the normal nuclear matter density 
(here, 0.16 fm~ 3 is used), p, p n and p p 
are the total, neutron and proton densi- 
ties, respectively. r 2 is zth component of the 
isospin degree of freedom, which equals 1 
or -1 for neutrons or protons, respectively. 
The coefficients a, (5 and 7 are parameters 
for nuclear equation of state (EOS). C sym 
is the symmetry energy strength due to the 
difference of neutron and proton. In the 
present work, we take a = -124 MeV, (3 = 
70.5 MeV and 7 = 2.0 which corresponds 
to the so-called hard EOS with an incom- 
pressibility of K — 380 MeV and C sym = 32 
MeV [25] . V c is the Coulomb potential and 
jjYuk j g Yukawa (surface) potential which 
has the following form: 



U Yuk = %!-Y— exp(Lm 2 ) 
2m tf. m V ; 
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with V y = 0.0074 GeV, m = 1.25 fm' 1 
and L = 2.16 fm 2 . The relative distance 
r ij = \ I- Experimental in- medium NN 
cross section parametrization which is en- 
ergy and isospin dependent is used in this 
work. 

The Pauli blocking effect in IDQMD model 
is treated separately for the neutron and 
the proton: whenever a collision occurs, we 
assume that each nucleon occupies a six- 
dimensional sphere with a volume of f? j2 
in the phase space (considering the spin 
degree of freedom), and then calculate the 
phase volume, V, of the scattered nucleons 
being occupied by the rest nucleons with 
the same isospin as that of the scattered 
ones. We then compare 2V/H 3 with a ran- 
dom number and decide whether the colli- 
sion is blocked or not. 



In the QMD model, the initial momentum 
of nucleons is generated by means of the 
local Fermi gas approximation. The local 
Fermi momentum is given by: 

F i F (?)=h^p i (f))\,(i = n, P ). (7) 

In the model, the radial density can be writ- 
ten as: 

I T 2 -\- T 2 L 

P( r ) = 2 7^77^72 ex P( ^F^): 
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The time evolution of the colliding sys- 
tem is given by the generalized variational 
principal. Since the QMD can naturally 
describe the fluctuation and correlation, 
we can study the nuclear clusters in the 
model [25,26,27,28,29,30]. In QMD model, 
nuclear clusters are usually recognized by 
a simple coalescence model: i.e. nucleons 
are considered to be part of a cluster if in 
the end at least one other nucleon is closer 
than r min < 3.5 fm in coordinate space and 
Pmin < 300 MeV/c in momentum space 
[25] . This mechanism has been extensively 
applied in transport theory for the cluster 
formation. 

Now we move to the calculations. For an 
example, we simulated 86 Kr + 124 Sn at 
25 MeV/nucleon and impact parameter of 
7-10 fm. 50,000 events have been accu- 
mulated. The systems tend to freeze-out 
around 120 fm/c. In this paper, we ex- 
tract the following physics results at 200 
fm/c. The upper panel of Fig. 1 shows 
transverse momentum dependence of ellip- 
tic flows for mid-rapidity light fragments. 
The range of transverse momentum for dif- 
ferent fragments is different according to 
their masses. From the figure, it shows that 
the elliptic flow is positive and it increases 
with the increasing p t . It reflects that the 
light clusters are preferentially emitted 
within the reaction plane, and particles 
with higher transverse momentum tend to 
be strongly emitted within in-plane , i.e. 
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stronger positive elliptic flow. In compari- 
son to the elliptic flow at RHIC energies, 
the apparent behavior of elliptic flow ver- 
sus p t looks similar, but the mechanism is 
obviously different. In intermediate energy 
domain, collective rotation is one of the 
main mechanisms to induce the positive el- 
liptic flow [31,32,33,34,35,36]. In this case, 
the elliptic flow is mainly driven by the 
attractive mean field. However, the strong 
pressure which is built in early initial ge- 
ometrical almond-type anisotropy due to 
the overlap zone between both colliding 
nuclei in coordinate space will rapidly 
transforms into the azimuthal anisotropy 
in momentum space at RHIC energies [11]. 
In other words, the elliptic flow is mainly 
driven by the stronger outward pressure. 
The lower panel displays the elliptic flow 
per nucleon function of transverse 

momentum per nucleon, and it looks that 
there exists the number of nucleon scaling 
when Pt/A < 0.25 GeV/c. This behavior is 
apparently similar to the number of con- 
stituent quarks scaling of elliptic flow ver- 
sus transverse momentum per constituent 
quark (pt/n) for mesons and baryons which 
was observed at RHIC [11]. 

The RHIC experimental data demon- 
strated a scaling relationship between 2- 
nd flow (v 2 ) and n-th flow (v n ), namely 

v n(Pt) ~ v 2 (Pt) [37]. It has been shown 
[38,39] that such scaling relation follows 
from a naive quark coalescence model [13] 
that only allows quarks with equal mo- 
mentum to form a hadron. Denoting the 
meson anisotropic flows by u n ,ju(Pt) and 
baryon anisotropic flows by v n> B(pt) , Kolb 
et al. found that v^uipt) = (l/4)u| Af(Pt) 
for mesons and v^sipt) = (1/3)^1 b(p*) 
for baryons if quarks have no higher-order 
anisotropic flows. Since mesons dominate 
the yield of charged particles in RHIC 
experimental data, the smaller scaling 
factor of 1/4 than the empirical value of 
about 1 indicates that higher-order quark 
anisotropic flows cannot be neglected. In- 
cluding the latter contribution, one can 
show that 
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Fig. 1. (a) Elliptic flow as a function of trans- 
verse momentum (pt)- Squares represent for 
neutrons, circles for protons, triangles for frag- 
ments of A = 2, diamonds for A = 3 and stars 
for A = 4; (b) Elliptic flow per nucleon as a 
function of transverse momentum per nucleon. 
The symbols are the same as (a). 
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where v n>q denotes the quark anisotropic 
flows. The meson anisotropic flows thus 
satisfy the scaling relations if the quark 
anisotropic flows also satisfy such rela- 
tions. However, this ratio is experimentally 
determined to be 1.2 [40], which means 
that the fourth-harmonic flow of quarks 
v\ must be greater than zero. One can 
go one step further and assume that the 
observed scaling of the hadronic t> 2 ac- 
tually results from a similar scaling oc- 
curring at the partonic level. In this case 
v\ = (f|) 2 and the hadronic ratio v^v\ 
then equals 1/4 + 1/2 = 3/4 for mesons and 
1/3 + 1/3 = 2/3 for baryons, respectively. 
Again, since this value is measured to be 
1.2, even the partonic v\ must be greater 
than simple scaling and quark coalescence 
models predict. 
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Recognizing the above behaviors of the 
flows at RHIC energies, we would like to 
know what the higher order momentum 
anisotropy in the intermediate energy is. So 
far, there is neither experimental data nor 
theoretical investigation for the higher or- 
der flow, such as t>4, in this energy domain. 
In the present work, we explore the behav- 
ior of t>4 in the model calculation. Fig. 2 
shows the feature of V4. Similar to the rela- 
tionship of V2/A versus pt/A, we plot V4/A 
as a function of Pt/A. The divergence of 
the different curves between different par- 
ticles in Fig. 2(a) indicates no simple scal- 
ing of nucleon number for 4-th momentum 
anisotropy. However, if we plot v^/A 2 ver- 
sus (pt/A) 2 , it looks that the points of dif- 
ferent particles nearly merge together and 
it means a certain of scaling holds between 
two variables. Due to a nearly constant 
value of V4/V2 in the studied p t range (see 
Fig. 2(c)) together with the number-of- 
nucleon scaling behavior of v-ijA vs pt/A, 
V4/A 2 should scale with (p t /A) 2 , as shown 
in Fig. 2(b). If we assume the scaling laws 
of mesons (Eq.(9)) and baryons (Eq.(lO)) 
are also valid for A = 2 and 3 nuclear clus- 
ters, respectively, then v^/v 2 for A = 2 and 
3 clusters indeed give the same value of 1 /2 
as nucleons, as shown in Fig. 2(c). Coin- 
cident ally the predicted value of the ratio 
of V4/V2 for hadrons is also 1/2 if the mat- 
ter produced in ultra-relativistic heavy ion 
collisions reaches to thermal equilibrium 
and its subsequent evolution follows the 
laws of ideal fluid dynamics [41]. It is inter- 
esting to note the same ratio was predicted 
in two different models at very different 
energies, which is of course worth to be 
further investigated in near future. Overall 
speaking, we learn that v^jv\ is approxi- 
mately 1/2 in nucleonic level coalescence 
mechanism, which is different from 3/4 for 
mesons or 2/3 for baryons in partonic level 
coalescence mechanism. 



Measurements of inclusive single-particle 
spectra from heavy-ion collisions indicate 
a simple empirical pattern of light frag- 
ment production: the observed invariant 
momentum-space density pa for fragments 
with mass number A closely follows the 
A-th. power of the observed proton density 
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Fig. 2. (a) v^/A as a function of pt/A for dif- 
ferent particles, namely, neutrons (squares), 
protons (circles), fragments of A = 2 (trian- 
gles), A = 3 (diamonds) and A = 4 (stars), (b) 
V4/A 2 as a function of (pt/A) 2 . (c) the ratios 
of V4/V2 for different particles vs pt- 

pf. It was shown that this power law be- 
havior is valid for spectra of participant 
fragments up to A = 14 with projectiles 
ranging from protons to Au at a variety 
of beam energies between 0.1 A and 15 A 
GeV [24]. However, there are rare exper- 
imental data and calculations to test the 
power law at intermediate energies. In the 
present work, we test the momentum-space 
power law up to A = 4 by IDQMD. 

The upper panels of Fig. 3 show the 
transverse momentum space densities 
p = A 2 dN/p t dp t . The results also show 
a level of adherence to power law behav- 
ior at high Pt/A, which is similar to what 
have been previously reported for single- 
particle-inclusive measurements [42], and 
reflects the persistence of momentum- 
space coalescence behavior for intermedi- 
ate energy 86 Kr + 124 Sn collisions. The 
middle panels of Fig. 3 depict the den- 
sity in the in-plane transverse momentum 
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Fig. 3. Upper panels: transverse momen- 
tum-space density for light fragments; Middle 
panels: in-plane transverse momentum space 
density; Bottom panels: azimuthal distribu- 
tions relative to the reaction plane. The open 
squares represent the density for protons, the 
solid circles for the fragments with A = 2, the 
open circles for the proton density squared 
normalized to the same area as the proton 
density, the solid triangles for ^4 = 3, the open 
triangles for protons to the power of 3, the 
solid diamonds for A = 4, the open diamonds 
for protons to the power of 4, and the stars 
for the A = 2 density squared. 

space (p x /A). The densities are also nor- 
malized to the same area as the proton 
density. Wang et al. first introduced the 
momentum-space power law to p(p x ) [42]. 
Again, our results show a level of adher- 
ence to power law behavior. The bottom 
panels of Fig. 3 depicts the azimuthal (0) 
distribution of fragments relative to the 
reaction plane. The factor k is chosen so 
that the mean values of dN/ d(f) are fixed to 
be 1. Essentially the power-law behavior 
remains for the azimuthal distributions of 
these light fragments. From all the above 
distributions, nucleon coalescence mecha- 
nism keeps valid in the momentum space 
densities. 



In summary, we applied IDQMD model 
to investigate the behavior of anisotropic 
flows, namely v<i and t>4, versus transverse 
momentum for the light fragments from 25 
MeV/nucleon 86 Kr + 124 Sn collisions at 7- 
10 fm of impact parameters. Both V2 and v,i 
generally show positive values and increase 
with p t . By the number-of-nucleon scaling, 
the curves of elliptic flow for different frag- 
ments approximately collapse on the simi- 
lar curve, which means that there exists an 
elliptic flow scaling on the nucleon number 
for light fragments. This phenomenon is 
similar to the NCQ scaling of elliptic flows 
of hadrons at RHIC energies where elliptic 
flow can be developed in the early partonic 
stage of collisions, but here it is a nucleonic 
level in intermediate energy collisions. For 
4-th momentum anisotropy V4, it seems to 
be scaled by v\ and v±/v\ ~ 0.5. It will 
be of very interesting if one can measure 
this ratio in intermediate energy HIC. In 
addition, we also investigated momentum- 
space densities of light fragments as func- 
tions of fragment transverse momentum 
p t , in-plane transverse momentum p x and 
the azimuth angle relative to the reac- 
tion plane. All these observables are well 
described by the momentum-space power 
law which has been well experimentally 
observed at high energies. All the above 
phenomena can be seen as an outcome 
of the nucleonic coalescence which results 
that the number-of-nucleon scaling of both 
flow and momentum space density in inter- 
mediate energy heavy ion collision. It may 
indicate that nucleonic matter may be one 
of the middle transit stage before chemical 
freeze-out takes place. 
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